Theory for Electron-Doped Cuprate Superconductors: 
rf-wave symmetry order parameter 
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Using as a model the Hubbard Hamiltonian we determine various basic properties of electron-doped 
cuprate superconductors like Ndg-ajCezCuCU and P^-ajCe^CuC^ for a spin-fluctuation-induced 
pairing mechanism. Most importantly we find a narrow range of superconductivity and like for hole- 
doped cuprates d x 2_ y 2 - symmetry for the superconducting order parameter. The superconducting 
transition temperatures T c (x) for various electron doping concentrations x are calculated to be much 
smaller than for hole-doped cuprates due to the different Fermi surface and a flat band well below 
the Fermi level. Lattice disorder may sensitively distort the symmetry d x 2_ y 2 via electron-phonon 
interaction. 
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One expects on general physical grounds if Cooper- 
pairing is controlled by antiferromagnetism that d-wave 
symmetry pairing should also occur for electron-doped 
cuprates 0. Until recently experiment did not 

clearly support this and reported mainly s-wave pairing 
p-pl. Maybe as a result of this, so far electron-doped 
cuprates received much less attention than hole-doped 
cuprates. Previously, we were rather successful in de- 
termining the doping dependence of antiferromagnetism 
in both electron- and hole-doped cuprates by using the 
Hubbard Hamiltonian ||. Applying this model to the 
hole-doped cuprates, many physical quantities like the 
normal-state pseudogap and the doping dependence of 
T c can also be described JTo| , pd| . 

Hence, to get an uniform theory we use here for the su- 
perconducting properties of electron-doped cuprates also 
as a model the 2D one-band Hubbard Hamiltonian 



H = - 



E 



cLc 



(1) 



Here, cj^ creates an electron with spin a on site i, U 
denotes the on-site Coulomb interaction, and ty is the 
hopping integral. For the optimally doped NCCO the 
dispersion efe and Fermi surface are taken in accordance 
with photoemission (ARPES) experiments (jl2|. Thus, 
we choose the parameters t = 138 meV and if — 0.30 in 
calculating 



-2t [cos k x + cos k y — 2t' cos k x cos k y + fi/2] 



(2) 



where the chemical potential [i describes the band filling. 
Here and in the following, we set the lattice constant 
a = b equal to unity. 

In Fig. ^ the results for e/. are shown. For comparison, 
the results of a tight-binding calculation with t = 250 
meV and if = 0, which is often used to describe the 
hole-doped superconductors, is also displayed. One im- 
mediately sees one important difference: in the case of 
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FIG. 1. Comparison of the energy dispersion tk for elec- 
tron-doped cuprates and for hole-doped ones. Data (open 
dots) are taken from Ref. The solid curve refers to our 

tight-binding calculation as described in the text. The dashed 
curve corresponds to Eq. (□) with t = 250 meV and if = 0. 



NCCO the flat band is approximately 300 meV below 
the Fermi level, whereas for the hole-doped case the flat 
band lies very close to it. Thus, one expects a smaller 
T c for electron-doped cuprates than for the hole-doped 
cuprates. Then, using in a spin-fluctuation- induced 
pairing theory in the framework of the so-called FLEX 
approximation |l^-|l5[], we calculate the doping depen- 
dence T c (x) and some other basic properties. 

In Fig. H we show results for the real part of the 
spin susceptibility at 100K in the weak-coupling limit 
for u> = (solid curve) and for u> — uj s f rs 0.47£ (dashed 
curve). u) s f denotes the spin fluctuation (paramagnon) 
energy, where a peak in Im x(Q? w ) occurs. The commen- 
surate structure of i?ex(q, lu = 0) is in accordance with 
recent calculations in Ref. [ fLsf , where it was pointed out 
that the exchange of spin fluctuations yield a good de- 
scription of the normal state Hall coefficient R# for both 
hole- and electron-doped cuprates. Furthermore, we also 
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FIG. 2. Momentum dependence of the real part of the spin 
susceptibility along the BZ route 

(0,0) -> (tt,0) -> (tt,7t) -> (0,0) at T = 100 K for uj = 
(solid curve) and w = u) s f ~ 0.47i (dashed curve). The main 
contributions to the corresponding pairing interaction come 
from q pa ir (along the anti-nodes) and Q pa ir (along the 'hot 
spots') as is illustrated in Fig. |H[ 



find a linear temperature dependence of the in-plane re- 
sistivity p a b(T), if we do not take into account an addi- 
tional electron-phonon coupling. This will be discussed 
later. Concerning the superconducting properties, it was 
stated in Ref. |l~7| that in contrast to the hole-doped su- 
perconductors the electron-doped systems may be also 
close to a d^-symmetry instability. However, within the 
picture of a spin-fluctuation-induced pairing this is defi- 
nitely not the case. Since the lower tiny peak favors d xy 
pairing symmetry and the dominating larger peak d x i_ y i 
symmetry (but is pair- breaking for dry-symmetry) , one 
understands why an underlying superconducting order 
parameter 0(k, uj) exhibits almost pure d x 2_ y 2 symme- 
try. 

In Fig. |j| we present our result for <fi(k, uj) for u> = 
and a doping x = 0.15 at T/T c — 0.8, where the gap 
has just opened. The gap function has clearly d x 2_ y 2- 
wave symmetry. This is in agreement with the reported 
linear dependence of the in-plane penetration depth for 
low temperatures and with phase-sensitive measure- 
ments j|. From our obtained result of a pure d x 2_ y 2- 
wave superconducting order parameter we expect a zero- 
bias conductance peak (ZBCP) |18| as observed for the 
hole-doped superconductors 0. However, its absence in 
the electron-doped cuprates may be attributed to small 
changes in the surface quality and roughness Jl9| ] or to 
disorder pof . Note, the incommensurate structure in the 
order parameter close to (tt, 0) results from the double 
peak structure in Re at w k uj s f = 0A7t shown in 
Fig. 0. Physically, it means that Cooper-pairing occurs 
not only for Q = (tt, tt), but mostly for uj = uj s j and for 
Q* = (tt — 5, tt + S). Furthermore, from Fig. 0, Fig. 0, 




FIG. 3. Calculated d x 2 _ !/ 2-wave symmetry of the super- 



conducting order parameter at T/T c 
first square of the Brillouin Zone. 



=0.8 for x=0.15 in the 



and Fig. ^|we conclude that no d X!/ -symmetry component 
is present in the superconducting order parameter, since 
the dominating d^^-type pairing suppresses d xy pair- 
ing. ARPES study might test this. 

In Fig. [| we present our results for the phase dia- 
gram T c (x) and Tn(x). We find that in comparison to 
hole-doped superconductors smaller T c values and super- 
conductivity in a narrower doping range as is also ob- 
served in experiment pl| ] . Responsible for this are poorer 
nesting properties of the Fermi surface and a flat band 
around (tt, 0) which lies well below the Fermi level. The 
narrow doping range is due to antiferromagnetism up to 
x = 0.13 and, for increasing x, rapidly decreasing nesting 
properties. We have calculated the Cooper-pair coher- 
ence length £o, i-e. the size of a Cooper-pair, and find 
similar values for electron-doped and hole-doped super- 
conductors (from 6 A to 9 A). If due to strong coupling 
lifetime effects the superfiuid density n s becomes small, 
the distance d between Cooper pairs increases. If for 
0.15 > x > 0.13 the Cooper-pairs do not overlap signif- 
icantly, i.e. d/£o > 1, then Cooper-pair phase fluctua- 
tions get important |23|,^4|,[ll]] . Thus we expect like for 
hole-doped superconductors T c oc n s . Assuming that n s 
increases approximately linearly from x~0.13 to x~0.15 
we estimate a T c which is smaller than calculated from 
</>(k, lo) = (see Fig. 4). As a consequence more exper- 
iments determining T c for x<0.15 should be performed 
to check on the Uemura scaling T c oc n s . 

The effect of electron-lattice coupling on superconduc- 
tivity should depend on lattice perturbations like oxy- 
gen deficiencies. Then, the isotope effect may show a 
distinct effect of electron-phonon coupling on T c . On 
general grounds we expect a weakening of the d x 2_ y 2- 
pairing symmetry if we include the electron-phonon in- 
teraction and if this plays a significant role. The ab- 
sence of an isotope effect (ao = d\n.T c /d\aM « 0.05) 
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FIG. 4. Phase diagram T(x) for electron-doped cuprates. 
The AF transition line is taken from Ref. ||. Inset: blow-up 
of the doping region 0.18 < x < 0.12. The solid curve corre- 
sponds to our calculated T c values obtained from </>(k, u) = 0. 
For a comparison, also experimental data are shown (squares 
from Ref. (22), circles from Ref. g, triangle from The 
dotted curve refers to T s tx n s . 



for doping x = 0.15 (see Ref. p7fl) suggests the presence 
of a pure d x 2_ y 2 -symmetry. We know from Fig. |^ that 
phonons connecting the Fermi surface with wave vector 
Qpair = C 71 "' 7r ) w ^ a dd destructively to the spin fluctua- 
tion pairing p8|] . If, due to exchange of spin fluctuations, 
a d x 2_ y 2 -symmetry instability is the dominant contribu- 
tion to the pairing interaction, an additional electron- 



phonon coupling with wave vector qp a j r = (0.57T, 0) 
would be also pair building. Note, we generally expect 
that due to the poorer nesting the pairing instability due 
to electron-phonon and spin fluctuation interaction be- 
come more easily comparable. In this case, the electron- 
phonon coupling would definitely favor s-wave symmetry 
of the underlying superconducting order parameter. This 
can be analyzed in detail by adding a term a 2 F(q,uj) to 
the pairing interaction [ p8[ . The corresponding phonon 
modes were calculated in Ref. [ p9| . Moreover, the inclu- 
sion of an electron-phonon interaction yields a quadratic 
term in the resistivity for lower temperatures ^q] as it is 
observed in experiment Q . 

To continue the discussion why the symmetry of the or- 
der parameter depends for electron-doped cuprates more 
sensitively on electron-phonon interaction, we show in 
Fig. H the calculated Fermi surface for optimally doped 
NCCO. Note, the topology of the Fermi surface for the 
electron-doped cuprates is very similar to optimally hole- 
doped Bi 2 Sr 2 CaCu208+5 (BI2212) as it was also pointed 
out recently in Ref. ]3l]j . We estimate that mainly 
no phonons are present along the edges (— 0.257T, n) — > 
(0.257T, it) bridging BZ areas, where the superconducting 
order parameter, $>(k, uS), is always positive (denoted by 
+/+). Note, attractive electron-phonon coupling bridg- 








FIG. 5. Calculated Fermi surface for (optimally doped) 
NCCO. The +(— ) sign and the dashed curve corresponds to 
the calculated momentum dependence (see Fig. ^) of the the 
d x 2_ y 2 gap function <f>(k,U) = 0) and its nodes, respectively. 



ing +/— areas (— 0,57r, — 0.57r) — ► (0.5-7T, 0.57r) is destruc- 
tive for d x 2_ y 2 Cooper pairing. Bowever, due to poorer 
nesting conditions, pairing transitions of the type +/— 
are somewhat contributing and then a mixed symmetry 
{d x 2_ y 2 + as} may occur. 

Further experimental study of the doping dependence 
of the oxygen-isotope effect are necessary for a better 
understanding of the role played by the electron-phonon 
interaction. For example, if due to structural distor- 
tion and oxygen deficiency in the CuCvplane the phonon 
spectrum F(q,u) changes significantly, then this affects 
ao and reduces T c . Possibly the reported large isotope 
effect of «o = 0.15 for slightly changed oxygen content, 
i.e. Ndi.85Ceo.i5Cu03.8, could be related to this p2| , ^3| . 
As an example, one might think of the oxygen out of 
plane B2« mode, which become active if O4 is replaced 
by O3.8 I!. 

In summary, our model for electron-doped cuprates 
yields like for hole-doped case pure d x 2_ y 2 symmetry 



pairing in a good agreement with recent experiments. 
In contrast to hole-doped superconductors, we find for 
electron-doped cuprates smaller T c values due to a flat 
band dispersion around (jr, 0) well below the Fermi level. 
Futhcrmore, superconductivity only occurs for a narrow 
doping range 0.18 > x > 0.13 because of the onset of an- 
tiferromagnetism, and, on the other side, due to poorer 
nesting conditions. We get 2A/fcsT c = 5.3 for x = 0.15 
in reasonable agreement with Ref. We argue that if 
the electron-phonon coupling becomes important, for ex- 
ample due to oxygen deficiency, then the s-wave pairing 
instability competes with c? a; 2_ y 2-wave symmetry. This 
might explain a possible s-wave order parameter as re- 
ported in earlier measurements. 
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